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ADAPTIVE  OPTICAL  FILTERING  TECHNIQUES 


The  purpose  of  this  study  Is  to  examine  tbe  potential  of 
using  optical  information  processing  technology  for  adaptive 
beamforming.  Tbe  adaptive  beamforming  problem  consists  of 
estimation  of  tbe  covariance  matrix  of  the  noise  field  and 
inversion  of  tbe  covariance  matrix  to  obtain  tbe  weights  for  tbe 
elements  of  tbe  antenna.  There  are  numerous  algorithms  that  have 
been  developed  for  this  problem.  Our  concern  here  is  to  examine 
tbe  adaptive  beamforming  problem  in  view  of  tbe  capabilities  of 
optics  and  identify  areas  where  optics  can  be  used  to  benefit. 

There  are  two  basic  approaches  to  the  problem:  iterative 
solutions  and  dlreot  oovariance  matrix  inversion.  In  the 
iterative  approacb,  tbe  oovarlaaoa  matrix  is  never  formed  and 
thus  a  direct  matrix  inversion  operation  is  not  required.  An 
optimum  set  of  weights  is  estimated  iteratively,  by  adjusting 
eaoh  weight  to  minimize  the  degree  of  correlation  between  the 
output  of  the  antenna  and  tbe  signals  received  at  the  individual 
elements.  The  implementation  of  the  algorithm  requires  an  array 
of  complex  multipliers  to  set  the  weights  and  an  array  of 
multipliers/accumulators  to  perform  tbe  required  correlation 
operation.  Simple  analog  optical  systems  can  be  configured  to 
implement  the  Iterative  algorithm.  For  instance  an  array  of 
multipliers  oan  be  implemented  by  an  array  of  LEDs  Imaged  through 
a  multichannel  acoustooptic  device.  Accumulation  of  analog 
products  oan  be  optically  realized  by  temporally  integrating  suoh 
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produots  on  a  temporally  integrating  deteotor  array  or  an 
optioally  addressed  spatial  light  modulator.  The  more  important 
issue  la  what  are  the  possible  advantages  of  the  optical 
Implementation  over  an  eleotronio  analog  solution  (the  approaoh 
that  is  primarily  used  nov).  The  aoouraoy  (which  is  direotly 
related  to  null  depth)  is  in  general  comparable  sinoe  both 
systems  are  analog  and  have  similar  limitations.  An  advantage  of 
the  optloal  system  is  smaller  size,  weight  and  possibly  oost. 
These  factors  may  be  quite  significant  in  a  space-borne 
application,  but  it  does  not  appear  to  be  significant  enough  to 
have  an  impact  for  airborne  or  ground  based  applications.  The 
difficulty  with  using  an  analog  optloal  system  to  perform 
adaptive  beamformlng  for  a  narrowband  linear  antenna  array  is 
that  the  optloal  system  is  asked  to  do  too  little;  only  a 
relatively  small  one  dimensional  array  of  multipliers  is  needed. 
Normally,  we  apply  optloal  signal  processing  technlquea 
effectively  to  problems  that  require  a  muob  heavier  processing 
load,  thereby  utilizing  the  full  oapability  of  optics  and  thus 
obtain  a  dear  advantage  over  alternate  approaches. 

There  are  two  possible  extensions  of  the  simple  one 
dimensional  adaptive  array  that  we  have  considered  thus  far,  to 
which  optloal  techniques  may  be  more  effectively  applied.  The 
first  is  large  two  dimensional  arrays  consisting  of  several 
thousand  elements  and  thus  requiring  an  equal  number  of  analog 
multipliers  for  adaptive  beamformlng.  In  this  case  the  ease  with 
whloh  suoh  large  parallel  arrays  of  analog  multipliers  oan  be 


optically  implemented  provides  «  clear  advantage.  A  second  area 
is  adaptive  beaaforalng  of  relatively  saaller  arrays,  but  for 
broadband  (spread  spectrum)  applications.  In  this  case,  effective 
interference  cancellation  can  be  achieved  only  by  having  multiple 
taps  at  each  antenna  element.  For  large  signal  bandwidths 
(several  hundred  MHz)  in  a  dense  Interference  environment  we  may 
require  that  each  antenna  element  have  hundreds  of  taps.  This 
again  leads  to  a  requirement  for  very  large  arrays  of  analog 
multipliers.  Adaptive  beamforming  of  broadband  antenna  arrays  is 
probably  the  most  promising  application  of  analog  optical  methods 
in  this  general  area.  We  have  studied  this  application  in  detail 
and  have  developed  optical  architectures  specifically  for  this 
problem.  A  detailed  description  of  this  processor  can  be  found  in 
the  papers  attached  as  Appendices  A  and  B  to  this  report.  Another 
approach  to  this  problem  is  described  in  RADC  report  RADC-TR-83- 
1561. 

The  most  obvious  limitation  of  the  analog  implementation  (be 
it  optical  or  eleotronlo)  is  limited  accuracy.  Moreover,  the  low 
acouraoy  restricts  the  algorithmic  flexibility  to  variations  of 
the  Iterative  algorithm,  which  has  the  added  problem  that  the 
convergence  time  is  a  function  of  the  eigenvalues  of  the  noise 
oovarlance  matrix.  Shown  in  Figure  1  is  a  diagram  reproduced 
from  RADC  report  RADC-TR-81-1 302  which  plots  the  null  depth  as  a 
funotlon  of  the  number  of  bits  of  accuracy  for  several  direct 
algorithms  as  well  as  the  iterative  algorithm.  It  is  seen  that 
for  the  Iterative  algorithm  an  increase  in  the  number  of  bits 


Cancellation  Ratio  (dB) 


beyond  eight  does  not  provide  significant  Improvement  in 
performance  beoauae  the  residual  error  due  to  the  finite  number 
of  iterations  beoomes  dominant.  Much  more  effective  interference 


suppression  can  be  obtained  with  a  direct  algorithm  if  sufficient 
acouraoy  is  used.  The  significance  of  this  result  when  we 
consider  methods  for  Improving  the  performance  of  an  optical 
adaptive  array  processor,  beyond  what  is  obtainable  with  an 
analog  system,  is  that  we  must  find  ways  to  drastically  improve 
the  acouraoy  and  also  find  ways  to  implement  direct  algorithms. 
It  does  not  appear  that  we  can  use  effectively  an  inaccurate 
optical  processor  that  implements  a  direct  algorithm  or  a  very 
aoourate  one  that  implements  an  iterative  algorithm.  A  large 
number  of  optical  architectures  have  been  developed  over  the  last 
several  years  for  performing  matrix  operations  (vector/matrix  or 
matrix/matrix  products).  These  processors  can  in  principle  be 
applied  to  the  implementation  of  direct  algorithms  by  structuring 
the  algorithm  so  that  it  can  be  implemented  by  successive  matrix- 
matrix  products.  Direct  algorithms  are  generally  algorithms, 

i.e.  the  number  of  arithmetic  operations  required  is  proportional 
to  the  oube  of  the  number  of  elements  in  the  array.  When  we 
implement  such  an  algorithm  with  successive  matrix/matrix 
produota,  normally  N  such  products  are  needed  and  since  eaoh 
matrix  operation  is  an  operation,  the  optical  implementation 
requires  operations  to  solve  the  problem.  Ve  have  examined 
the  possibility  of  optioally  implementing  a  matrix  inversion  with 
]|3  operations  but  we  were  unable  to  identify  a  viable  solution. 
If  such  an  implementation  oan  be  identified,  then  this  would  be  a 


gnlfioant  advance  since  it  would  become  much  easier  to  obtain  a 


•eed  advantage  with  optics.  We  have  also  examined  the  problem 
1  pipelining  successive  matrix  operations  as  is  required  for  the 
iplementatlon  of  a  direct  matrix  algorithm.  A  space  integrating 


itrix-matrix  multiplier  was  developed  that  has  the  capability  of 
illy  efficient  pipelining:  the  format  of  the  product  matrix  is 
:actly  the  same  as  that  of  one  of  the  input  matrices.  This 


rstem  is  described  in  the  paper  attached  as  Appendix  C  to  the 
iport. 

The  accuracy  of  essentially  all  analog  matrix  processors 
in  be  improved  by  processing  digitized  (binary)  data  through  the 
»e  of  the  digital  multiplication  by  analog  convolution  algorithm 
iMAC).  This  algorithm,  can  provide  accuracy  superior  to  that  of 
ie  analog  implementation,  but  it  reduces  the  processing  power 
;he  maximum  number  of  mu 1 t 1 p 1 1 ca t ions/addi t ions  that  can  be 
iplemented)  and  Increases  the  complexity  of  the  post  detection 


Lectronlcs,  since  analog  to  digital  converters  are  required  to 
iduce  the  signal  produced  by  the  optical  system  to  true  binary 
>rmat.  In  what  follows  we  examine  these  trade-offs. 

A  schematic  diagram  of  a  generalized  optical  processor  that 
irforms  linear  operations  using  the  DMAC  algorithm  is  shown  in 
lg.2.  The  system  has  input  parallel  spatial  channels  and  each 


tie  of  these  channels  accepts  binary  bits  at  a  rate  B  y .  Each 
bannel  may  accept  bits  from  a  separate  external  information 
>urce  or  alternatively,  each  channel  accepts  the  bit  that  was  at 


tie  adjacent  channel  during  the  previous  cycle  (one  oyclezl/B^ 


V\.v--vv\ 


soonds).  The  information  in  each  channel  is  multiplied  by  M 
Lfferent  bits  in  the  optical  processor.  We  call  M  the  fan-out 
actor.  The  system  has  N2  parallel  output  channels  each  having 


amporal  bandwidth  B2.  The  signal  detected  at  each  output  channel 
s  then  electronically  converted  to  the  binary  representation  by 
n  A/D  converter.  Having  defined  these  parameters  we  can 
haracterize  the  performance  of  any  specific  architecture  without 
arther  knowledge  of  the  details  of  the  implementation.  Thus  we 
ill  be  able  to  derive  some  guidelines  that  are  generally 
pplicable.  The  number  of  ■ b  i  t  -  mu  1 1 i p 1 i c a t i o n a"  that  the 
rooessor  performs  per  unit  time  is  equal  to  N^B^M.  The  number  of 
it  multiplications  that  are  required  to  realize  one 
ul tiplication  between  two  Integers  with  DMAC  is  aK2,  where  "a* 
s  a  oonstant  between  1  and  U,  depending  on  the  efficiency  of  the 
peciflc  Implementation,  and  K  is  the  number  of  bits  that  are 
sed  to  represent  eaoh  number.  The  processing  power,  P,  of  the 
verall  system  is: 

P  s  N^B^M/aK2  multiplications  per  second.  (1) 

1  ear 1 y ,  P  is  a  number  we  wish  to  maximize.  The  number  of  analog 
o  digital  conversions  that  need  to  be  performed  per  unit  time 

a : 

C  s  N2B2  A/D  conversions  per  second.  (2) 

gain  it  is  clear  that  C  is  a  number  we  wish  to  minimize  in  order 
.o  keep  the  complexity  of  the  electronicsat  a  minimum. 
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jure  2.  Generalized  Optical  Matrix  Processor  Utilizing  the 
DMAC  Algorithm. 


ortunately,  P  and  C  cannot  be  independently  chosen.  To 
reciate  this  fact,  we  define  the  ratio 


R  =  P  /  C  multiplications  per  A/D  conversion.  (3) 

ncreases  monotonical ly  as  either  P  increases  or  C  decreases, 
refore  we  want  to  make  R  as  large  as  possible.  R  provides  us 
h  a  direct  comparative  estimate  of  the  optical  implementation 
'sus  an  all  digital  implementation.  If  for  instance  R  were 
al  to  one,  then  only  one  binary  multiplication  is  being 
formed  by  the  optical  system  per  A/D  conversion.  Since  it  is 
>ut  equally  difficult  to  perform  multiplications  and  A/D 
iverslons  electronically,  this  would  be  a  strong  indication 
it  optics  offers  no  advantage  in  this  case.  In  order  to 
ermine  the  maximum  value  for  R  we  consider  the  characteristics 
the  output  stage  of  the  processor.  The  number  of  bits  that  are 
Lng  generated  per  unit  time  within  the  optical  processor  is 
lal  to  aK2P.  The  number  of  samples  that  are  being  transferred 
:  of  the  optical  processor  per  unit  time  is  C.  The  ratio 
(2P/C)  is  equal  to  the  maximum  number  of  bits  that  is 
iumulated  (through  either  spatial  or  temporal  integration)  to 
•m  each  output  sample.  Therefore,  this  ratio  cannot  exceed  the 
;put  dynamic  range  of  the  system,  DR2,  which  is  defined  as  the 
jber  of  distinguishable  signal  levels  that  are  produced  by  the 
tector  and  A/D  converter.  This  gives  us  a  very  simple  upper 
ind  for  R: 

R  <  DR2/aK2. 
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(4) 


of  the  array  and  also  la  the  temporal  reference  signal  used  for  the  detection  of 
ed  signal,  s(t).  It  can  t>e  derived  by  feeding  a  tapped  delay  line  with  e(t);  the 
om  the  ith  tap  corresponds  to  si( t>  and  the  tap  spacing  is  adjusted  to  obtain  the 
>ok  direction: 

■  s(t-ia)  «  s(t-i(d/c)cosa) ,  (J) 

s  the  desired  look  angle  from  boresight,  d  is  the  array  element  spacing,  and  c  is 
i  of  light. 


Fig. *  Multi-Channel  Acoustooptic  Correlator 

\\\ 

1  1  k 


proc«j»or  OUTPUT 


Fig. 10  Broadband  Adaptive  Array  Processor 

shows  the  optically  implemented  adaptive  array  processor  with  the  AOD-impl emented 
me  filter  in  the  upper  branch  and  the  n-channel  correlator  in  the  lower  one.  By 
I  Eqs.(7)  and  <*),  the  equation  that  determines  the  filter  function,  h^ft)  is  seer 


wrwrrf?»?*r 


i  c  a  1  processors  are  particularly  attractive  candidates  for  adaptive  broadband  array 
1  processing  because  of  their  inherent,  multidimensional  processing  capability, 
irer.  the  broadband  requirement  can  be  met  by  the  use  of  broadband  multi-channel  AOD's; 
arge  time-bandwidth  products  that  are  available  translate  to  the  possibility  of 
nenting  very  high  order  transversal  filters.  The  acoustooptic  processor  which  we  now 
Lbe  is  an  extension  of  the  active  temporal  processor  described  earlier,  to  the  space- 
domain.  Here,  we  employ  a  combination  of  multi  and  single  channel  AOD's  to  perform 
squired  operations. 


e  output  of  a  general,  space-time  filter,  with  a  finite  accumulation  time  can  be 

ised  as 


5  T/2 

t>  -  >  /  dx, 

i-1 


(4) 


xi(t)  is  the  signal  from  the  ith  antenna  element,  and  h^(t>  is  the  filtering  function 
ch  channel.  A  similar  expression  is  obtained  for  the  optically  implemented  space-time 
r  using  two  multi-channel  AOD's  shown  in  Fig.*.  The  output  of  this  filter  is 


t) 


N  t/2 

2  ( 1 /T) /  xi(t-T)h1 (t+x)  dt, 
i-l 


-T/2 


(5) 


*e  see  that  the  only  difference  from  the  general  filter  (Eq.4)  is  in  the  time 
ession  of  the  output.  It  can  be  shown  that  the  optimum  choice  for  h^(t)  in  Eq.(S) 
fies  the  system  of  integral  equations: 


T/2 


(«) 


/  h* (t  +  X ’ ) T < X-X ' )  dt '  -  XSj(t-X), 

-T/2  J  J 

!  8j(t)  is  the  desired  signal  vector  and  T^jtt)  is  the  covariance  matrix  of  the  input 
[  given  by  J 


j(x)  -  E  fni(t)nj(t-x ) J . 

is  the  noise  vector  appearing  at  the  array  elements. 


(7> 


>r  adaptivity,  we  need  to  calculate  and  continuously  update  the  filter  function  h^(t) 
rive  the  output  to  the  optimum  result.  As  with  the  previously  described  Active 
>E80r.  the  output  must  be  correlated  with  the  input  to  produce  the  appropriate  filter 
tion.  Since  the  array  processor  has  n  inputs  and  one  output,  this  requires  that  we 
slate  n  signals  with  a  common  one.  This  can  be  achieved  with  the  arrangement  shown  in 
1  which  shows  the  use  of  a  multi-channel  AOD  driven  by  the  n  antenna  element  outputs, 
onjunction  with  a  single  channel  AOD  which  is  driven  by  the  array  output  signal, 
tfically.  the  output  of  the  ith  element  of  the  linear  detector  array  of  the  correlator 
.ven  by 


T/2 


i(t>  -  (i/T)  /  y < t*x ) x j ( t»2x ) dx . 
-T/2 


<*) 


proper  correlation  to  appear  at  each  output,  the  signal  driving  the  single  channel  AOP 
be  time-compressed  by  a  factor  of  two.  This  is  indeed  the  case  for  the  system 
ribed.  and  thus,  the  AOD  implemented  space-time  filter  and  the  n-channel  correlator 
a  single  reference  are  compatible. 


hown  in  Fig.io  is  the  array  processor  system  diagram  that  shows  the  interconnections 
are  required;  it  is  a  direct  extension  of  the  Active  Processor  to  2-dimensions.  The 
ut  from  each  antenna  element  is  correlated  with  the  processor  output  to  produce  the 
er  function  for  that  element.  The  steering  vector,  s^(-t).  determines  the  look 
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Pig.*  Adaptive  Optical  Procasaor  (Active) 

Illj.  Adoptive  ACQUStfiCEtiC  AlIOV  £l££i£SCX 

general  adaptive  acheme  is  shown  in  Fig. 7  in  which  each  channel  corresponding  to  one 

>  element  is  filtered  by  a  transversal  filter  with  n  controllable  tap  weights;  if  the 

>  has  N  antenna  elements,  then  the  array  processor  requires  the  adaptive  control  of  nK 
tendent  weights.  Therefore,  the  computational  load  for  broadband  adaptive  phased  array 
easing  is  much  heavier  than  that  for  adaptive  temporal  filtering  because  here,  the 
r  processor  must  adapt  in  two  independent  dimensions:  space  and  time. 
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Correlation  can  be  performed  with  the  convolution  configuration  just  described  by  time- 
everaing  the  input  to  one  of  the  AOD's.  This,  however,  is  not  acceptable  for  real-time 
peration.  Shown  in  Pig. 4  is  a  apace-integrating  processor  which  provides  the  relative 
lotion  between  the  two  input  aignala  that  is  required  for  correlation  by  using  a 
lemagnifying  system  between  the  two  AOD's.  The  diffracted  light  from  the  first  AOD  is 
•aged  with  a  2:1  magnification  onto  the  second  AOD.  and  the  diffracted  light  from  both 
>OD‘s  is  spatially  integrated  onto  a  single  detector.  If  the  input  to  the  second  AOD  is 
iae-compressed  by  a  factor  of  two,  i.e..  i4(t)  ■  i$(2t).  then  the  correlator  yields 

T/4  t+T/2 

I«(t)  *  I  ij(t*2*)i4(t+t)dT  “  /  i.<T)ij(t  +  OdT.  (J) 
-T/4  t-T/2 

is  its  output.  The  above  is  seen  to  be  a  finite  window  correlation  of  the  two  signals. 
j(t)  and  i 5 ( t > .  The  advantage  of  obtaining  the  time  compressed  output  from  the  convolver 
s  now  clear,  since  it  can  be  used  as  the  input  to  the  correlator  just  described  and  thus 
ibtain  a  consistent,  cascaded  convolution-correlation  operation. 


The  optical  implementations  of  the  Passive  and  Active  Processors  are  shown  in  Figs.  S 
ind  4.  respectively,  and  they  differ  only  in  their  electrical  interconnections.  Fig.  5(6) 
Is  the  system  of  Fig. 1(2)  with  the  optically  implemented  blocks  of  Figs.  3  and  4.  In  both 
Eases,  the  upper  branch  of  the  processor  computes  the  correlation  while  the  lower  performs 
the  convolution.  Since  the.  convolver  and  correlator  have  a  common  input,  the  first  AOD  is 
ihared.  A  more  detailed  description  can  be  found  in  the  previous  publication4. 
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Fig.i  Active  Processor 


For  signal  detection,  where  a  known  broadband  signal  is  to  be  detected  amidst  an 
idditive  collection  of  narrowband  jammers  of  unknown  frequencies,  a  different  system  is 
required.  Shown  in  Pig.  1  is  a  system  diagram  of  such  a  system,  the  Active  Processor,  along 
sith  a  brief,  mathematical  description  of  its  operation;  a  complete  description  can  be 
found  in  reference  4.  As  with  the  Passive  Processor,  the  basic  operation  of  this  system  can 
»e  understood  using  heuristic  arguments.  Suppose  that  the  input  consists  soley  of  a  large 
implitude  sinusoid  of  frequency  f«.  and  the  desired  waveform.  s(t).  is  a  broadband  signal, 
rith  the  loop  opened  at  the  summing  junction,  the  output.  e(t).  will  consist  mostly  of  a 
linusoid  of  frequency  f«;  this  signal  is  then  correlated  with  the  input  to  produce  another 
linusoid  of  the  same  frequency  at  the  summing  junction.  When  the  loop  is  closed,  the 
>utout  sinusoid  resulting  from  the  convolution  of  the  input  and  the  code  waveform.  s(t).  is 
iso  out  of  phase  with  that  which  is  produced  by  the  feedback  signal,  and  hence, 
rancellation  of  the  sinusoid  takes  place.  Again,  the  feedback  signal  is  negligible  for 
imall  amplitude  broadband  signals. 

Both  the  Passive  and  Active  systems  require  convolution  and  correlation  as  the  basic 
>uilding  blocks.  In  choosing  the  architecture  to  be  used,  we  consider  the  following 
Important  requirements:  Dwide  bandwidth.  J )dynamical ly  controllable  convolver  and 
sorrelator.  Dconvolver  -  correlator  compatibility.  4>large  dynamic  range.  These 
considerations  lead  us  to  choose  a  space-integrating  acoustooptic  architecture. 
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rig.l  Space-Integrating  Acoustooptic  Convolver 

Convolution  can  be  performed  with  two  AOD's  as  shown  in  Fig. 3.  The  diffracted  signal 
from  the  first  AOD  is  imaged  with  unity  magnification  onto  the  second  and,  in  the  coherent 
realization  shown,  the  diffracted  components  are  spatially  integrated  onto  a  single 
detector  whose  output  photocurrent  is 

T/I 

I.(t)  ■  /  i.ft+Oijft— Odr.  t«*/v,  T»W/2.  <1) 

-T/a 

l|(t)  and  i2<t)  are  the  complex  amplitudes  of  the  AOD  input  signals,  v  is  the  acoustic 
velocity  in  the  AOD.  and  V  is  the  physical  length  of  the  AOD  aperture.  This  integral  can  be 
•anipulated  to  yield  the  more  familiar  form: 

t*T/2 

I,(t)  •  /  1](y)1n( at  — t ) dr ,  (2) 

t-T/a 

which  is  recognised  to  be  a  finite  window  convolution,  time-compressed  by  a  factor  of  two. 
The  seemingly  troublesome  compression  is  actually  an  advantage  in  this  case  as  will  be 
apparent  when  we  consider  the  correlator  implementation. 
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Space-integrating,  acoustooptic  processors  for  adaptive,  temporal  filtering  are 
examined.  The  basic  architecture  is  then  extended  to  the  space-tiae  domain  for  application 
in  broadband  phased  array  processing.  An  acoustooptic  processor  capable  of  such  2- 
diaensional.  adaptive  processing  is  described. 
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A  aajor  portion  of  optimum  filtering  theory  concerns  itself  with  the  efficient 
separation  of  useful  signals  from  additive  noise.  Fixed  optimum  filters,  such  as  the 
wiener  filter,  are  applicable  when  the  signal  and  noise  statistics  are  stationary  and  known 
apriori;  the  lack  of  such  apriori  knowledge,  however,  motivates  the  implementation  of 
adaptive  filters  which  estimate  the  required  signal  and  noise  characteristics.  The 
implementation  of  such  filtering  techniques  requires  a  processor  which  must  be  able  to 
compute  various  correlation  functions  of  signal  and  noise  and  change  its  transfer  function 
accordingly. 

The  situation  becomes  considerably  more  complex  when  one  is  required  to  adaptively 
filter  signals  in  the  space-time  domain,  as  in  the  case  of  broadband  array  processing.  The 
number  of  broadband  jammers  that  an  adaptive  processor  can  cancel  without  compromising  its 
directional  discrimination  can  be  used  as  its  performance  measure.  Electronic 
implementations  of  such  array  processors1  have  exhibited  limited  performance  with  respect 
to  this  measure.  The  transversal  filters  responsible  for  temporally  filtering  the  outputs 
from  the  array  elements  are  of  low  order  and  hence,  operation  is  usually  limited  to 
narrowband  signals,  due  to  hardware  limitations. 

The  bandwidth  requirements  and  the  parallel  nature  of  array  processing  make  optical 
Implementations  attractive.  Various  optical  implementations  have  been  explored  by 
researchers  in  the  area2'*.  In  a  recent  paper4,  we  described  two  optical  adaptive  filters 
for  use  in  the  time  domain.  These  implementations  are  strictly  one-dimensional.  requiring 
only  one  dimensional  devices,  permitting  an  extension  to  the  space-time  domain  through  the 
use  of  multi-channel  AOD's.  In  this  paper,  we  briefly  review  the  operation  of  the 
temporally  adaptive  acoustooptic  filter.  After  this  foundation  has  been  established,  an 
adaptive  array  processor  utilizing  multi-channel  AOD’s  in  a  space  integrating  architecture 
will  be  described. 


II*  AdABliXt  ZAB&QXAl  CiitAXA 

To  be  adaptive,  a  processor  that  optimizes  either  the  mean-square  error  or  SNR  criterion 
must  compute  various  correlation  functions  of  the  input  signal  and  noise  and  vary  its 
filter  response  characteristics  accordingly.  Specifically,  we  will  consider  situations 
where  the  signal  is  broadband  and  the  additive  noise  consists  of  strong,  narrowband  jammers 
whose  frequencies  are  unknown;  the  jammer  frequencies  must  be  estimated  for  effective  noise 
rejection. 

Shown  in  Fig.)  is  a  system  diagram  of  the  Passive  Processor  which  was  shown  to 
adaptively  perform  an  approximate  Wiener  filtering  operation4.  The  operation  of  the  systen 
can  be  explained  heuristically  in  the  following  manner.  Suppose  that  the  input  consists 
soley  of  a  single  sinusoid,  and  consider  the  feedback  signal.  x(t),  with  the  loop  opened  at 
the  summing  junction.  The  input  sinusoid  is  first  convolved  with  a  replica  of  itself, 
delayed  by  •  to  produce  another  sinusoid  of  the  same  frequency,  also  delayed  by  «.  at  the 
input  of  the  correlator.  This  is  then  correlated  with  the  input  delayed  by  •  to  produce  a 
sinusoid  with  the  same  phase  as  the  input  sinusoid  appearing  at  the  summing  junction;  the 
delays.  •.  cancel  out  due  to  the  cascaded  correlation-convolution  operations.  The 
possibility  of  stable  cancellation  is  now  apparent  since  the  two  signals,  the  input  and 
feedback  signals,  which  are  subtracted  at  the  junction,  are  identical  in  frequency  and 
phase.  Cancellation  does  not  occur  for  a  broadband  input,  because  it  correlates  poorly  with 
delayed  versions  of  itself,  resulting  in  a  negligible  feedback  signal.  fi(t).  The  passive 
processor  thus  discriminates  against  narrowband  components  of  the  input  signal,  while 
preserving  the  desired,  broadband  components  of  the  input,  making  it  suitable  for  signal 
estimation. 
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procaasor  to  ibtiaguish  batwsmi  tha  various  jamming 
atgnnh  passant  at  the  input 
Tbs  optksDy  implsmsntad  actlv  procnssof  la  thown 
achaaaaticaMy  in  Fig.  8;  tha  active  system  diffata  from 
tha  pmaiva  only  in  ha  electrical  interconnections.  Hart 
tha  output  of  tha  ayatam  ia  a  turn  compressed  signal 
since  it  ia  taken  at  tha  output  of  tha  canvoKsr,  this 
presents  no  problems  tinea  tha  output  signal  is  a  cor* 
relation  peak,  and  we  are  inUraatad  in  detection  rather 
than  estimation.  As  in  tha  passive  ease,  tha  convolver 
output  drives  AOD3  of  the  correlator.  Tha  correlator 
output  ia  than  subtracted  from  the  time  reversed  ref¬ 
erence  signal  (the  matched  filter  impulse  response)  to 
form  tha  input  to  AOD2  of  tha  convolver,  note  that  even 
though  tha  output  ia  time  compressed,  the  time  scaling 
ia  compatible  for  all  tha  signals  within  the  system.  To 
see  this  more  precisely,  as  examine  tha  input-output 
aquation  for  the  active  system;  using  Eq*.  (16)  and  (17), 
we  obtain 

e(l)  -  Um  f™iU  +  t)o»(r-tHr 
J-TH 

-  (G/T*)  «(t  ♦  -  r  ♦  Se)r(t  -  r  ♦  eMedr. 

<»> 

where  r  (t )  is  the  system  output.  Under  low  SNR  con¬ 
ditions,  Eq.  (23)  can  be  approximated  by 

/•r/4 

«(()  a  (1/D  I  i(l  +  »)«*(f-lWT 
J-TH 

-  (C/T*)  n(f  ♦  r  )•»•(!  -  t  ♦  S»)e(i  -  f  ♦  v)4udr, 

(14) 

The  output  spectrum  is  easily  calculated  from  Eq.  (24) 
to  be 

(%e)J*  X(n)S*(w  -  •)  Hi*  ”  Sn)T/4vJe 
IM  e  — — — — — — — — — — —  • 

1  ♦  (C/4 a*)  J*’  |N(«)|*Usc*|(m  -  U)T/*o\da 

m 

Aside  from  the  appearance  of  the  sine  functions  which 
arise  due  to  the  finite  apertures  of  the  AODs  and  the 
time  compression,  Eq.  (25)  is  quite  similar  to  Eq.  (9),  the 
matched  filter  result 

V.  Conclusion 

Two  adaptive  systems  particularly  suitable  for  use 
in  narrowband  jamming  noise  environments  have  been 
presented.  The  passive  system  has  been  shown  to  be 
appropriate  for  estimation  purposes  while  the  active, 
for  signal  detection.  A  particularly  interesting  exten¬ 
sion  of  this  work  is  the  implementation  of  time  -space 
adaptive  filters  for  processing  broadband  signals  from 
phased  array  antennas.  The  processors  dweribed  ia 
this  paper  are  temporal  filters,  but  since  they  are  im¬ 
plemented  in  only  one  dimension,  H  is  possible,  through 
use  of  multiple  transducer  AODs,  to  implement  a  2-D 
time  space  adaptive  filter.  The  1-D  processors  are 


currently  being  built  in  the  laboratory,  and  future  re¬ 
search  wiU  include  the  experimental  verification  of  the 
operation  of  the  systems  as  weU  as  statistical  charac¬ 
terization  of  tha  effects  of  the  finite  apertures  of  the 
AODs  and  detector  noise. 
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2(u)  a  XM  -  (G/4**)  JJ’  N‘(a) 

X  N(fi)  npl~j{0  -  a)a]Z(u  +  a  -  fi) 

X  «inc|(w  -  «)77 4»]  Hnc|(w  +  a  -  25)T/4r)darf0.  (20) 

The  double  integral  of  (20)  can  be  approximated  by  a 
•ingle  integral  since  it  becomes  appreciably  large  only 
near  the  region  a  »  0,  so  that 

ZM  »  X(«>  -  <C/4**)Z(w)  J**  |N(a>|» 

X  mdc*](w  -  a)T/4w\da.  (21) 

Solving  for  Z(w)  gives 

MOO  muB>  OPTICS  /  Vet ».  Me.  tf  /  1 0e«cH  1SS4 


v.vv.v.  vv  -.vV  v. 


1  +  (C/4**)  J]]  |N(a)|* tinc*((w  -  a)774*|do 

(22) 

Except  for  the  presence  of  the  sine  function  in  the  de¬ 
nominator  due  to  the  finite  apertures  of  the  AODs,  Eq. 
(22)  is  quite  similar  to  the  Wiener  result,  Eq.  (6). 
Shown  on  Figs.  6  and  7  are  plots  of  Z{n>)tX\w)  tot  an 
input  consisting  of  two  sinusoids,  one  at  9.8  MHx  and 
the  other  at  10.2  MHz;  the  first  plot  is  for  T  ■  10  ^sec, 
while  the  second  is  for  T  ■  50  msoc.  Clearly,  the  finite 
aperture  of  the  AODs  limit  the  resolving  power  of  the 
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been  incorporated  into  one  system  requiring  only  one 
light  source  and  three  AODs.  The  correlator  and  con¬ 
volver  share,  in  effect,  the  AOD  which  is  driven  by  the 
system  input  signal  x  ( t ).  The  input  signal,  translated 
in  frequency  too^is  applied  to  AODl  and  the  summing 
node  before  entering  AOD2.  The  convolution  is  per¬ 
formed  by  the  lower  leg  of  the  processor  and  the  corre¬ 
lation  by  the  upper  leg.  As  shown  in  the  previous  sec¬ 
tion,  the  optical  convolver  takes  its  two  inputs,  each  at 
the  carrier  frequency  u>o  and  gives  their  convolution, 
compressed  in  time  by  a  factor  of  2  and  translated  to 
2u>o,  «  output.  This  output  is  then  used  to  drive 
AOD3  of  the  correlator,  recalling  that  to  achieve  a 
proper  correlation,  one  of  the  correlator  inputs  must  be 
time  compressed  by  a  factor  of  2  and  translated  to  2«o, 
we  see  that  its  output  is  free  of  time  scaling  problems. 
This  signal  is  then  subtracted  from  the  input  signal  to 
form  the  system  output  z(t)  which  is  fed  back  to  the 
convolver  via  AOD2.  The  delay  element  which  is  nec¬ 
essary  in  the  passive  processor  can  be  realized  simply 
by  translation  of  AODl  along  the  direction  in  which  the 
acoustic  signal  propagates. 

We  now  derive  the  output  of  the  optical  passive 
processor  that  is  obtained  under  the  conditions  stated 
in  Sec.  II.  We  begin  by  stating  the  mathematical  op¬ 
erations  performed  by  the  optically  implemented  blocks 
of  Pigs.  1  and  2;  with  <i  and  ij  as  the  inputs  of  the  blocks, 
the  convolver  block  output  is  given  by 

Jrn 

Htf  ♦  -  rWr.  (16) 

■Tl* 

and  the  correlator  block  output  is  given  by 


i*(t)  -  (l IT)  C™  ij(t  +  2t )Mt  +  r)dr.  (17) 

*/-TY4 

Equations  (16)  and  (17)  are  simply  Eqs.  (13)  and  (14), 
respectively,  with  the  carrier  frequencies  suppressed. 
These  integals  are  different  from  their  counterparts 
found  in  Eqs.  (1)  and  (2)  in  two  important  ways.  First, 
Eqs.  (16)  and  (17)  have  finite  limits  of  integration  lim¬ 
ited  by  the  maximum  transit  delay  time  of  the  AODs. 
Also  the  arguments  of  the  integrand  are  different,  and 
these  differences  reflect  the  time  compression  problems 
of  the  optical  convolver  and  correlator  referred  to  ear¬ 
lier.  In  both  cases,  the  integration  limits  extend  from 
-774  to  Tl 4,  where  T  is  the  aperture  of  each  AOD  in 
units  of  time.  With  Eqs.  (16)  and  (17),  the  input-out¬ 
put  equation  of  the  passive  processor  is 

lit)  -  x(t)  -  (C/T»)  Jj  *•(*  +  2v-v) 

X  xit  +  r  +  v  -  *)*(t  -  r  +  vMrrfi v,  (18) 

where,  as  in  Sec.  II,  z(f )  is  the  system  output,  and  x{t) 
-*(()  +  n(t)  is  the  input;  »(f)  is  the  broadband  signal, 
and  n(t)  is  the  narrowband  interference.  The  delay 
element  a  causes  the  broadband  components  to  become 
decorrelated  so  that  the  narrowband  components 
dominate  the  double  integral  of  Eq.  (18)  yielding 

i(t)  •*(!>-  (C/T*)  n»(t +  -  ») 

X  nit  ♦  r  ♦  v  —  »)*(!  —  r  ♦  vUlrdv  (19) 

Fourier  transformation  then  yields 

1  October  1M4  /  Vd.  23.  Mb.  19  /  APPUED  OPTICS  3479 


A- 5 


of  the  AOD;  noU  that  the  Bragg  condition.  Eq.  (10),  waa 
used  to  derive  the  laat  line  of  Eq.  (11).  The  firat  term 
on  the  right-hand  side  of  Eq.  (11)  corresponds  to  the  dc 
or  undiffracted  light  which  is  unmodulated  to  first 
order.  The  second  term  corresponds  to  the  (-1)  dif¬ 
fracted  order,  which  is  modulated  by  aj(t  -  x/v),  Dop- 
pler  shifted  by  -«o,  and  deflected  by  an  angle  20b-  The 
amplitude  distribution,  Ai(x,t)  in  Eq.  (11),  is  then  im¬ 
aged  onto  the  second  AOD  by  lenses  L 1  and  L2\  the 
imaging  reverses  the  spatial  coordinate  z.  The  second 
AOD  is  positioned  at  the  Bragg  angle  with  respect  to  the 
dc  component  from  AODl,  and,  therefore,  a  portion  of 
it  is  diffracted  at  an  angle  20b  by  AOD2.  The  major 
portion  of  the  diffracted  beam  from  AODl  passes 
through  AOD2  undiffracted,  since  it  is  not  Bragg- 
matched,  and  it  propagates  in  the  same  direction  as  the 
beam  diffracted  by  AOD2.  The  remaining  components 
of  the  light  at  the  exit  of  AOD2  are  angularly  separated 
from  these  two  components  and,  therefore,  can  be 
blocked  at  the  focal  plane  of  lens  L3.  The  two  compo¬ 
nents  of  interest  are  given  by 

Atlxjt)  «  [s2<<  -  x/v)  np{-j2r  liafatA)  tspOW) 

+  ai<t  +  x/v)  exp <-;'2v  sinSpxA) 

X  expt-juoOl  rect(x/tV).  (12) 

The  first  term  on  the  right-hand  side  of  Eq.  (12)  is  the 
(+1)  diffracted  order  from  AOD2,  which  is  modulated 
by  a2<t  -  x/v),  Doppler  shifted  by  wo,  and  deflected 
from  the  dc  by  20b;  this  beam  is  collinear  with  the  (-1) 
diffracted  beam  from  AODl  [Le.,  the  second  term  in  Eq. 
(12)].  The  light  amplitude  at  the  back  focal  plane  of 
lens  L3  is  the  Fourier  transform  of  Eq.  (12).  A  detector 
with  an  active  area  sufficiently  large  to  integrate  the 
entire  transform  is  placed  at  the  Fourier  plane,  and  the 
resulting  photocurrent  is 

/»<*>  -  |  _  |S|(I  +  r)|*dr  +  f  |a,U-r)|»dr 
•f-Tn  %f-rn 


Fig.  4.  Spacs  integrating  optical  comlator. 

coordinate  of  the  signal  in  AOD3,  so  that  the  two  signals 
are  copropagating  at  AOD4.  The  imaging  system  has 
a  2:1  demagnification  ratio,  so  that  the  velocity  of  the 
image  of  the  light  that  is  diffracted  by  AOD3  is  half  of 
that  of  the  acoustic  signal  in  AOD4.  Thus  the  signals 
in  the  AODs  continuously  translate  with  respect  to  each 
other.  The  undiffracted  beam  from  AOD3  is  incident 
at  the  Bragg  angle  of  AOD4  corresponding  to  its  2wo 
center  frequency.  The  demagnification  by  a  factor  of 
2  increases  the  angular  separation  between  the  dif¬ 
fracted  and  undiffracted  beams  from  AOD3  at  the  plane 
of  AOD4.  Therefore,  the  two  diffracted  beams  from 
AOD3  and  AOD4  are  parallel  to  one  another  when  they 
enter  lens  L6.  16  forms  the  Fourier  transform  of  the 
light  exiting  AOD4,  and  a  photodetector  spatially  in¬ 
tegrates  the  entire  spectrum.  Through  a  development 
similar  to  that  illustrated  for  the  convolver,  we  find  that 
the  photocurrent  from  the  detector  in  Fig.  4  is 
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+  J  Rc  np(/W)  J_n  »i«  ♦  r)t,(l  -  rWf  ,  (13)  +  2  Rc  eipOW)  J  +  2r)a,<r  ♦  rWr  .  (14) 


where  r  ■  x/v,  and  T  is  the  acoustic  delay  through  the 
aperture  of  the  AOD,  i.e.,  T  “  W/v.  The  first  two  terms 
on  the  right-hand  side  of  Eq.  (13)  are  low  frequency 
components,  and  they  are  removed  by  filtering.  The 
third  term  is  the  desired  operation  which  is  the  convo¬ 
lution  of  aid)  and  a2(t).  The  convolution  is  com¬ 
pressed  in  time  by  a  factor  of  2  and  translated  to  twice 
the  original  carrier  freqeuncy.  Notice  that  in  both 
passive  and  active  processors,  the  time  compressed 
output  of  the  convolver  is  one  of  the  inputs  to  the  cor¬ 
relator.  Thus  the  correlator  must  be  designed  to  accept 
a  time  compressed  signal  (by  a  factor  of  2)  in  one  of  its 
input  ports. 

In  contrast  to  convolution  which  requires  an  inversion 
of  a  coordinate,  correlation  requires  that  one  of  the 
input  signals  be  in  motion  with  respect  to  the  other  with 
no  coordinate  inversion.  Referring  to  Fig.  4,  the  cor¬ 
relator  consists  of  two  identical  AODs  with  counter- 
propagating  acoustic  signals;  u>o  end  2u>o  are  the  carrier 
frequencies  of  the  signals  driving  AOD3  and  AOD4, 
respectively.  The  imaging  system  reverses  the  spatial 
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where  k  is  the  spatial  frequency  variable  at  the  output 
plane.  The  limits  of  integration  of  the  correlator  are 
half  of  those  of  the  convolver.  This  implies  that  the 
apertures  of  the  AODs  used  in  the  correlator  must  be 
twice  as  long  as  those  in  the  convolver  to  obtain  equal 
integration  limits.  Unilluminated  portions  of  the  AODs 
in  the  correlator  and  the  convolver  can  be  used  for  in¬ 
troducing  delays.  The  third  term  is  the  desired  oper¬ 
ation  which  resembles  a  correlation  translated  to  the 
carrier  frequency  w o.  If  t2<t)  *  aa(2£),  the  third  term 
of  Eq.  (14)  becomes 

eipOW)  (  +  a)da|  ,  (1$) 

Jl-Tt*  \ 

which  is  a  finite  aperture  correlation  of  the  signals  ai(t) 
and  «3(t)  placed  on  the  carrier  frequency  wp. 

IV.  Adaptive  Optical  Processor  Implementations 

A  schematic  diagram  of  the  optically  implemented 
passive  processor  is  shown  in  Fig.  5.  Through  use  of  a 
beam  splitter,  both  the  convolver  and  correlator  have 
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The  Fourier  transform  of  the  output  signal  e(t)  is 
readily  found  to  be 
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1  ♦  C|JV{«)|*  C|iVW>J* 
where  £M  is  the  Fourier  transform  of  e(t),  and  the  last 
result  is  valid  if  the  feedback  gain  G  is  made  sufficiently 
large.  If  the  input  SNR  is  high,  it  is  evident  from  the 
equations  that  the  system  will  equalize  the  spectrum  of 
the  total  input  signal  while  simultaneously  performing 
a  matched  filtering  operation.  Identifying  |N(w)|2  as 
the  power  spectrum  of  the  noise,  Eq.  (9)  is  the  matched 
filter  result.20  Thus,  from  Eqs.  (6)  and  (9),  both  active 
and  passive  processors  share  the  common  property  of 
suppressing  strong  narrowband  components  that  are 
present  in  the  input  signal.  This  property,  commonly 
called  line  cancellation,  is  suitable  in  environments 
where  sinusoidal  jammers  are  present  for  either  esti¬ 
mation  or  detection  of  broadband  spread  spectrum 
signals. 


M.  Optical  Filters 

The  two  adaptive  filters  that  we  described  in  the 
previous  section  are  both  implemented  with  a  convolver 
and  a  correlator.  There  are  several  possible  ways  to 
implement  an  optical  correlator/convolver.  In  selecting 
an  implementation  that  is  best  suited  for  this  applica¬ 
tion,  we  must  consider  the  following  requirements:  (a) 
The  optical  filters  should  be  capable  of  processing 
broadband  signals  (1-GHz  bandwidth)  to  be  applicable 
to  spread  spectrum  systems.  This  requirement 
suggests  an  acoustooptic  implementation,  (b)  The 
impulse  response  of  the  convolver  and  the  correlator 
must  be  dynamically  controllable.  The  impulse  re¬ 
sponse  of  the  correlator,  in  particular,  is  also  a  broad¬ 
band  signal,  (c)  The  two  filters  must  be  compatible 
with  one  another;  the  impulse  response  of  the  convolver 
is  determined  by  the  output  of  the  correlator,  and  the 
output  of  the  convolver  is  one  of  the  inputs  to  the  cor¬ 
relator.  (d)  We  must  select  architectures  that  have  the 
maximum  linear  dynamic  range  possible.  Adaptive 
filters  are  typically  used  to  process  signals  with  very  low 
SNRs;  this  implies  that  they  must  have  a  sufficient 
dynamic  range  to  place  very  deep  nulls  at  the  frequen¬ 
cies  where  the  interference  occurs  in  order  to  suppress 
it  effectively.  The  acoustooptic  space  integrating 
convolver,  consisting  of  two  counterpropagating  AODs, 
an  integrating  lens,  and  a  single  detector,  is  well  suited 
for  this  application  in  that  its  two  inputs  and  the  output 
are  broadband  electrical  signals,  and  a  very  high  dy¬ 
namic  range  is  possible  in  this  architecture.  One  po¬ 
tential  difficulty  with  using  this  type  of  convolver  is  the 
fact  that  its  output  is  time-compressed  by  a  factor  of  2, 
and  thus  it  cannot  be  directly  correlated  with  the  input 
signal  as  required  for  adaptive  filtering.  A  second 
problem  is  that  correlation  cannot  be  readily  performed 
with  this  architecture;  one  of  the  input  signals  must  be 
time  reversed  before  it  is  applied  to  the  convolver  to 
compute  the  correlation  with  such  a  system.  These  two 
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Fig.  3.  Space  integrating  optical  convolver. 


problems  combine  to  provide  a  solution  in  this  case. 
We  will  show  that  a  space  integrating  correlator  can  be 
implemented  with  two  copropagating  AODs  if  one  of  the 
inputs  to  the  correlator  is  time  compressed.  The  output 
of  the  space  integrating  convolver  can  thus  be  used  as 
one  of  the  inputs  to  such  a  correlator.  In  this  section  we 
describe  a  specific  coherent  realization  of  the  acous- 
tooptic  convolver  and  correlator  blocks.  Several  vari¬ 
ations  of  these  processors  are  possible,  including  an 
incoherent  implementation.  We  begin  with  the  de¬ 
scription  of  the  convolver  shown  in  Fig.  3.  It  consists 
of  two  AODs  with  copropagating  acoustic  signals.  The 
1:1  imaging  system  between  the  AODs  provides  the 
coordinate  inversion  of  one  signal  with  respect  to  the 
other  which  is  necessary  for  the  convolution.  Let  s  i  it ) 
and  *2 (t)  which  are  of  the  form,  sit)  *  ait)  cos(u>o*  + 
pit)}  "  Re(a(7)  exp(/u>oO],  represent  the  signals  applied 
to  the  AODs;  the  complex  envelope,  ait)  -  ait) 
expljpit )],  is  the  modulating  signal,  and  wo  is  the  center 
frequency  of  the  AODs.  The  first  AOD  in  Fig.  3  is  il¬ 
luminated  by  a  collimated  coherent  light  beam  incident 
at  the  Bragg  angle  #b  (Ref-  21): 


»in0B  “  Xu>o/4vu, 


(10) 


where  X  is  the  wavelength  of  light  in  the  acoustooptic 
crystal,  and  u  is  the  speed  of  sound  in  the  crystal. 

The  interaction  between  the  acoustic  wave  induced 
by  the  signal  SiU)  and  the  optical  field  causes  some  of 
the  incident  light  to  be  diffracted,  and  since  the  AODs 
are  operated  in  the  Bragg  regime,  only  first -order  dif¬ 
fraction  is  appreciable.  For  a  more  compact  analysis, 
the  induced  acoustic  signal  in  the  AOD  is  treated  as  the 
system  input;  i.e.,  the  fixed  delay  between  the  applied 
electrical  signal  and  the  acoustic  signal  is  ignored.  The 
final  results,  however,  are  strictly  equivalent.  In  the 
actual  implementation,  additional  delays  will  be  nec¬ 
essary  to  compensate  for  the  delays  intrinsic  to  the  ac¬ 
oustooptic  system.  For  weak  modulation,  the  optical 
field  at  the  exit  plane  of  AODl  is  given  bv22 

/t(x.f)  *  |np(-;2v  »in0*i/A) 

♦  ;lm/2>s',(t  -x/c)  -*A')J 

x  eip(-;2v  tinltpx/XH  rvcUx/lV) 

«  |expl-;2*  sin9»x/X)  ♦ /(ro/2lal(t  -  x/il  »ip<-;uiof) 
x  **p(j2»  *in((sx/X)|  rectlx/H'),  (11) 

where  m  is  a  constant,  x  is  along  the  direction  of  the 
acoustic  wave  propagation,  and  IV  is  the  aperture  size 
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presence  of  additive  noise,  but  it  requires  that  the  au¬ 
tocorrelations  as  well  as  the  cross-correlation  of  the 
signal  and  the  noise  be  known  in  advance.  The  filter 
that  maximizes  the  output  SNR  is  the  matched  fill  r 
which  requires  prior  knowledge  of  the  autocorrelation 
of  the  noise  and  also  the  desired  signal  waveform. 

When  the  necessary  correlation  functions  are  not 
known  a  priori,  they  must  be  estimated  adaptively  from 
the  past  history  of  the  signal  and  the  limited  and  qual¬ 
itative  information  that  is  available.  We  assume  that 
the  available  information  is  as  follows:  the  signal  and 
noise  are  uncorrelated,  and  the  signal  is  broadband,  as 
in  spread  spectrum  systems.  We  model  the  wideband 
signal  as  bandiimited  white  noise,  so  that  S,(w),  the 
spectral  density  of  the  wideband  signal  s(t),  is  equal  to 
So  rect((«>  -  wo)/Aw],  where  wo  is  the  center  frequency 
of  the  signal  and  Aw  is  its  bandwidth.  The  spectral 
density  of  the  interference  n(t)  is  assumed  to  be  un¬ 
known,  and  it  is  adaptively  estimated. 

Shown  in  Fig.  1  is  a  system  diagram  of  the  passive 
processor,  which  will  be  shown  to  adaptively  perform 
an  operation  that  approximates  the  Wiener  filter  under 
the  conditions  stated  in  the  above  paragraph.  The  two 
necessary  operations  in  implementation  of  this  filter  are 
a  convolution  and  a  correlation.  The  convolver  serves 
as  the  variable  filter  that  is  controlled  by  the  correlator 
which  estimates  the  correlation  of  the  interference. 
The  operation  of  the  convolver  and  correlator  Mocks  are 
described  by  the  following  input-output  relations: 


u(t) »  J*  i|(l  -  rXitrWr. 

«) 

u.-d)  ■  J*  »»<!  +  rli'jdWr, 

(2) 

where  *,«)  and  ij(t)  are  the  inputs  to  the  blocks.  Using 
Eqs.  (11  and  (2),  we  find  the  relationship  between  the 
input  and  output  signals,  z(t)  and  z(t),  respectively,  of 
the  passive  processor  of  Fig.  1: 

*(<>  •  j rtf)  —  C  Jpf*  j(r)r*tar  -  c)xit  ♦  a  -  r  -  a)dadr 

■  xd)  -  C  ^  x<r)|*»<a  -  o)  ♦  n*(a  -  »)| 

x  |«<(  ♦  a  -  r  -  a)  +  nil  +  or  -  t  -  n)]dadr,  (3) 

where  a  is  a  time  delay.  Low  input  SNR,  along  with  the 
assumption  of  the  uncorrelatedness  between  the  signal 
tit )  and  the  noise  n(t),  results  in  the  reduction  of  Eq. 
(3)  into 

xd)  x(r)n*(a  -  *)nd  +  a  -  r  -  o)dadi. 

(4) 

Fourier  transformation  of  the  above  equation  yields  the 
approximate  steady-state  description: 

ZM  a  XM  -  GZMI/VM*.  (5) 

where  Z(w)  and  X  («)  are  the  Fourier  transforms  of  r  (t ) 
and  x(t),  respectively,  and  |N(w)|2  is  the  power  spec¬ 
trum  of  a  sample  realization  of  n(t ).  Solving  for  Z(w) 
gives 
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The  input  signal  to  the  convolver  and  the  correlator  in 
Fig.  1  is  delayed  by  the  same  amount,  and,  therefore,  the 
output  of  the  convolver-correlator  combination  i(t) 
becomes  independent  of  the  delay.  This  is  an  impor¬ 
tant  result  since  it  assures  that  the  feedback  signal  f  (f ) 
is  inphase  with  the  input  signal,  and  hence  the  system 
will  be  stable.  With  G  “  1/S0and  identifying  |N(w)|2 
as  the  estimate  of  S„  (w),  the  spectral  density  of  n  (t ),  Eq. 
(6),  describes  the  output  of  a  Wiener  filter.20 

We  now  show  that  the  active  processor,  shown  in  Fig. 
2,  adaptively  performs  matched  filtering.  The  exact 
input-output  relationship  is  given  by 


. /: 


xd  -  r)»*(-rWr 


“Cjj*  xd  -  r)x*(«)4(r  +  atdadr,  (7) 

where  e(t )  is  the  output  signal  of  the  active  processor. 
The  assumption  of  low  input  SNR  along  with  the  con¬ 
ditions  stated  previously  allows  the  following  approxi¬ 
mate  form  of  Eq.  (7): 


c(i)  a  J*  xd  -  r UV-rldr 

-  C  Jfj*  nit  -  t)«»(«)4(t  ♦  a\dadr 


(8) 
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Adaptive  acoustooptic  fitter 

Appendix  A 

Demetri  Psaltis  and  John  Hong 


A  new  adaptive  filter  utilitinc  acoustooplic  device*  in  a  (pace  intefratinc  architecture  i*  deacribed.  Two 
configuration*  are  presented,  one  of  them,  auitable  for  eigne)  estimation,  is  shown  to  approximate  the  Wie¬ 
ner  filter,  while  the  other,  suitable  for  detection,  is  shown  to  approximate  the  matched  filter. 


L  Introduction 

The  design  of  optimum  systems,  in  the  classical  sense, 
requires  a  priori  some  knowledge  of  the  signals  to  be 
encountered.  As  a  result,  such  systems  perform  poorly 
when  the  appropriate  characteristics  of  the  input  signals 
are  not  known  a  priori  sufficiently  well  or  are  time- 
varying.  An  adaptive  processor  has  the  ability  to  self- 
optimize  by  continually  monitoring  its  performance  and 
updating  its  parameters.  Adaptive  techniques  have 
been  applied  to  both  spatial  and  temporal  filtering  do¬ 
mains.  Specifically,  adaptive  techniques  have  been 
applied  to  antenna  array  processing  by  Appelbaum,1 
Widrow,2  and  others.3-4  Applications  to  time  domain 
problems  include  Lucky’s  work  on  data  redundancy 
removal,3  Sondhi’s  adaptive  echo  canceler,6  Widrow's 
work  on  adaptive  noise  suppression,7  Morgan  and 
Craig's  adaptive  linear  predictor,8  and  more.8-11 

In  the  work  on  time  domain  problems  cited  above  the 
adaptive  filtering  scheme  is  based  on  the  orthogonality 
principle.12  The  basic  idea  behind  the  scheme  is  to 
control  a  variable  filter  so  as  to  minimize  the  correlation 
between  the  input  signal  and  the  residual  signal,  which 
is  the  difference  between  the  input  signal  and  the  filter 
output.  A  particular  implementation  of  this  scheme 
uses  the  transversal  filter  architecture,  which  consists 
of  a  tapped  delay  line,  variable  filter  weights,  and  a 
summer  that  produces  a  weighted  sum  of  delayed  ver¬ 
sions  of  the  input  signal  as  its  output. 

The  linearity  and  parallel  nature  of  the  transversal 
filter  arrangement,  commonly  called  the  correlation 
cancellatiuirioop  (CCL)  system,  make  the  optical  im¬ 
plementation  possible.  The  advantages  of  optical 
processors  in  terms  of  bandwidth  and  the  large  effective 
number  of  taps  make  such  an  implementation  attrac¬ 
tive.  For  example,  acoustooptic  devices  (AODs),  which 
can  serve  as  optically  tapped  delay  lines,  are  superior 
in  terms  of  bandwidth  to  CCD  implemented  delay  lines 
that  are  currently  used  in  analog  adaptive  filters.1314 
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Several  adaptive  optical  filter  implementations  have 
been  previously  proposed.  Psaltis  et  al. 1516  proposed 
use  of  an  iterative  electrooptic  processor  for  adaptive 
spatial  filtering  of  phased  array  antenna  signals. 
Rhodes17  described  a  system  using  acoustooptic  and 
electrooptic  modulators  in  a  time  integrating  architec¬ 
ture  to  implement  the  CCL  algorithm  in  the  time  do¬ 
main,  and  VanderLugt18  described  an  optical  processor 
which  is  a  frequency  domain  implementation  of  the 
CCL  algorithm.  Lee  et  al. 19  have  devised  an  adaptive 
filter  which  suppresses  narrowband  interference  from 
wideband  signals  using  an  acoustooptic  spectrum  ana¬ 
lyzer  with  an  array  of  electrooptic  modulators  in  the 
spatial  frequency  plane  which  can  adaptively  excise 
strong  narrowband  components  of  the  signal  spec¬ 
trum. 

In  this  paper,  we  describe  a  new  method  for  adaptive 
inverse  filtering  in  the  time  domain.  Two  implemen¬ 
tations  are  presented;  the  first,  the  passive  processor, 
performs  the  inverse  filtering  operation,  while  the  other, 
the  active  processor,  is  capable  of  performing  adaptive 
matched  filtering  operations.  Both  processors  are  space 
integrating  systems  and  do  not  require  optically  ad¬ 
dressed,  time  integrating  spatial  light  modulators. 
These  architectures  implement  the  CCL  algorithm  in 
the  time  domain  and  are  structurally  simpler  than  their 
frequency  domain  implemented  counterparts.  The 
implementations  are  strictly  one  dimensional,  and  thus 
it  is  possible  to  extend  this  concept’ to  2-D  processing 
problems  through  the  use  of  multichannel  AODs. 

The  adaptive  processors  are  discussed  at  the  system 
level  in  Sec.  II,  and  the  optical  implementations  of  the 
system  blocks  are  described  in  Sec.  III.  The  two  optical 
configurations  of  the  adaptive  filters  are  described  in 
Sec.  IV. 

II.  Adaptive  Estimation  and  Detection 

The  extraction  of  information  from  a  signal  corrupted 
by  additive  noise  requires  a  priori  knowledge  of  the 
properties  of  the  desired  signal  and  the  noise.  If  the 
necessary  information  is  known,  the  optimum  linear 
filters  can  be  designed  to  satisfy  the  specified  perfor¬ 
mance  criteria.  The  Wiener  filter  provides  the  mini¬ 
mum  mean  square  error  estimate  of  a  signal  in  the 
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For  example  If  K=10  and  a=1,  we  require  the  number  of 
distinguishable  levels  at  the  output  to  be  at  least  100  for  R  to 
be  equal  to  1.  It  is  important  to  note,  that  the  output  levels 
must  be  sufficiently  well  defined  so  that  the  A/D  converter  can 
detect  all  of  them  with  very  low  probability  of  error.  DR2  is 
much  smaller  than  what  is  conventionally  called  "detector  dynamic 
range”.  We  estimate  that  it  would  require  very  sophisticated  - 
engineering  to  obtain  DR2=100  and  it  does  not  appear  that 
DR2s1000  is  practically  feasible  in  the  foreseeable  future.  We 
believe  that  the  answer  to  obtaining  a  high  accuracy  optical 
processor  that  is  suitable  for  implementation  of  direct 
algorithms  for  adaptive  antenna  arrays  is  not  likely  to  come  from 
a  well  engineered  D MAC  based  system.  It  is  more  likely  that 
algorithms  and  architectures  can  be  found  that  utilize  more 
efficiently  the  capabilities  of  optics  and  rely  less  on  the  use 
of  sophisticated  post  detection  electronics. 

1  RADC-TR-83-1 56,  "Acoustooptic  Adaptive  Processing  (AOAP)",  by 
General  Electric  Co.,  Dec.  1983. 

2  RADC-TR-8  1  -  1  30,  "Efficient  Digital  Algorithms  for  Adaptive 
Arrays",  by  Niel  A.  Carlson,  1981. 
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For  comparison.  consider  the  Fourier  Transform  of  Eg.*;  when  the  accumulation  time  is 
infinite,  the  optimum  filter  equation  is  given  by 

N 

2  J  f  ^  j  (•)  -  A*sJ(v),  (12> 

j-1 

Tijie)  •  F.T.  C  J. 

where  f t -i <•* >  is  the  spectral  density  matrix,  identifying  the  integral  in  Eq.il  as  the 
smoothed-’eatlmate  of  the  spectral  density  matrix  of  the  input  noise  vector.  Eq.  11  is 
approximately  equivalent  to  Eq.12.  The  effect  of  the  finite  time  integration  window  is 
seen  in  the  smoothing  of  the  noise  spectrum. 


v: 


♦  tin;  vector 


x(t)  ■  output 


M-.r  -K-Mm: 


Fig. 11  Adaptive  Acoustooptic  Phased  Array  Processor 
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An  optica!  architecture  is  described  for  performing  pipelined  matrix-matrix  multiplications.  The  architecture  is  imple¬ 
mented  using  multiple  transducer  acousto-optic  devices  and  a  wideband  photodetector  array.  A  variant  of  engagement  for¬ 
matting  allows  multiple  inner  products  to  be  simultaneously  computed  by  1  -D  spatial  integration,  and  through  proper  pipe¬ 
lining  the  full  product  matrix  is  produced  at  the  output  of  the  detector  array.  The  output  matrix  in  this  architecture  is  in  a 
format  that  is  directly  compatible  with  the  input ,  a  feature  that  can  facilitate  the  implementation  of  iterative  matrix  algo¬ 
rithms.  Digital  multiplication  by  analog  convolution  can  be  incorporated  for  improved  accuracy  by  using  a  frequency  multi¬ 
plexed  representation  of  the  binary  dau. 


1 .  Introduction 


The  optical  implementation  of  matrix  operations 
has  received  considerable  attention  recently.  Architec¬ 
tures  and  algorithms  have  been  designed  that  have  in¬ 
creased  the  speed,  accuracy  and  flexibility  of  optical 
matrix  processors,  extending  the  potential  applicabili¬ 
ty  of  such  systems  to  a  broader  range  of  problems. 
Specific  advances  that  have  been  accomplished  in  re¬ 
cent  years  include  the  initial  demonstration  of  vec  or- 
matrix  multiplication  [1 J,  the  introduction  of  time  in¬ 
tegrating  systolic  [2], engagement  [3,4], and  outer 
product  optical  processors  [5], a  frequency  multi¬ 
plexed  processor  [6], improvements  in  accuracy  with 
residue  arithmetic  [7],  the  utilization  of  the  method 
of  digital  multiplication  by  analog  convolution  (DMAC) 
(8,9]  in  the  above  array  processors  [10,1 1  ],  and  a 
combination  of  systolic  processing  and  the  DMAC  al¬ 
gorithm  in  a  two  dimensional  implementation  utilizing 
crossed  multichannel  Brag  cells  for  matrix  vector  mul¬ 
tiplications  with  digital  accuracy  [12]. 

Perhaps  the  most  significant  application  of  numeri¬ 
cal  optical  processors  is  in  (XAj)  problems,  i.e.  ma¬ 
trix  algebra  problems  that  require  a  minimum  of  M3 
multiplications  and  additions,  where  A  is  the  size  of 
the  matrix  involved.  The  solution  of  a  set  of  linear 
equations,  matrix  inversion,  and  singular  value  decom¬ 


position  are  examples  of  such  problems  [1 3].  Optical 
techniques  can  be  applied  to  such  problems  by  select¬ 
ing  an  algorithm  that  can  implement  the  required 
operation  with  successive  matrix-matrix  multiplica¬ 
tions,  such  as  Gauss  eliminations,  Givens  transforma¬ 
tions  or  Housholder  reflections  [13].  For  these  algo¬ 
rithms  A'  optical  matrix  multiplications  are  required, 
and  since  each  matrix-matrix  multiplication  requires 
A3  multiplies  and  adds,  optical  systems  usually  solve 
an  A 3  problem  with  A4  operations.  However,  the 
speed  and  parallelism  of  optics  can  make  the  optical 
implementation  advantageous,  despite  this  inefficiency. 
The  product  matrix  that  is  produced  at  each  iteration 
during  the  execution  of  such  an  algorithm  is  used  as 
one  of  the  input  matrices  for  the  next  iteration.  It  is 
therefore  important  that  the  formal  of  the  output 
product  matrix  is  directly  compatible  with  the  input 
in  order  to  avoid  reformatting  and  minimize  the  itera¬ 
tion  time.  The  architecture  described  in  this  paper  was 
selected  principally  because  the  output  can  be  ampli¬ 
fied  and  applied  directly  to  the  input.  A  space  inte¬ 
grating  implementation  using  a  parellel  output  wide- 
band  photodetector  array  is  chosen  for  accuracy  and 
speed  considerations.  Several  candidate  data  flow  opti¬ 
cal  architectures  satisfying  these  requirements  are  pos¬ 
sible.  In  this  paper  we  present  one  such  data  flow  ma¬ 
trix  processor  which  uses  crossed  multichannel  acousto- 
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optic  devices  (AOD).  The  operation  of  the  system  is 
based  on  time  and  space  alignment  of  vectors  which 
allows  the  formation  of  multiple  inner  product  sum¬ 
mations  via  spatial  integration  in  a  pipelined  fashion. 
The  matrix  formal  is  similar  to  an  engagement  array, 
but  the  data  flow  is  transposed  sq  that  the  local  mul¬ 
tiplications  needed  for  each  inner  product  operation 
form  in  parallel  in  space,  rather  than  sequentially  in 
time.  Each  inner  product  is  summed  by  a  1-D  space 
integrating  condensing  cylinder  onto  an  output  detec¬ 
tor.  Many  such  inner  product  accumulators  are  multi¬ 
plexed  in  the  orthogonal  dimension  onto  separate  de¬ 
tectors  of  a  linear  array. 


2.  Optical  processor  architecture 

The  proposed  optical  architecture  for  matrix-ma¬ 
trix  multiplication  is  shown  in  fig.  1 ,  along  with  the 
appropriate  data  flow.  The  principal  components  of 
the  system  are  a  pulsed  laser,  two  orthogonal  multi¬ 
channel  Bragg  cells  (one  with  N  channels,  and  the 
other  with  2A'  -  1  channels),  a  linear  array  of  A'  wide¬ 
band  photodetectors,  and  lenses.  The  optical  processor 
is  a  2-D  array  of  A'2  analog  multipliers  configured  as  an 
array  of  N  space  integrating  inner  product  processors. 
With  the  appropriate  engagement  format  of  the  ma¬ 


trices  an  array  of  inner  products  is  formed  on  the  de¬ 
tector  array  during  each  processor  cycle.  As  data 
flows  through  the  Bragg  cells  the  output  appears  in  the 
same  engagement  format  as  the  input,  which  allows  di¬ 
rect  feedback  for  iterative  operations  without  latency. 

Global  system  synchronization  and  sample  defini¬ 
tion  at  the  output  are  provided  at  each  time  interval 
T,  by  the  strobing  action  of  a  repetitively  pulsed  laser 
diode .  The  pulsed  light  is  collimated  and  incident  on 
the  first  multichannel  acousto-optic  device ,  AOD1 ,  at 
the  Bragg  angle  in  the  x  dimension.  The  elements  of 
an  A'  X  N  matrix  A  are  applied  to  the  A'  transducers 
of  AOD1  in  an  engagement  representation.  They  prop¬ 
agate  continuously  along  the  x  direction  at  a  velocity 
equal  to  one  inter  transducer  spacing  of  AOD2  each  T 
s.  Rows  are  represented  in  individual  channels  as  se¬ 
quential  acoustic  pulses  separated  by  T  s.  The  nth  row 
of  the  matrix  A  is  applied  to  the  nth  transducer  of 
AOD1  with  a  delay  of  nT  s.  In  this  manner  the  matrix 
is  folded  back  in  time  into  a  sliding  parallelogram  for¬ 
mat  we  call  time  engagement.  The  optical  field  emerg¬ 
ing  from  AOD1  is  spatially  filtered  in  the  Fourier 
plane  to  remove  the  undiffracted  component  and  the 
diffracted  field  is  imaged  onto  AOD2  at  the  Bragg 
angle  iny>.  For  clarity  the  image  reversal  of  the  imag¬ 
ing  system  is  ignored.  Matrix  B  propagates  in  AOD2 
in  the  y  direction  one  channel  separation  of  AOD1 
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Fig  J.  Schematic  representation  of  the  space  integrating  acousto-optic  matrix-matrix  multipber.  with  the  associated  engagement 
data  flow .  (For  simplicity ,  the  details  of  spatial  filtering  and  the  effect  of  image  reversal  arc  neglected  tn  this  figure  > 
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each  ft.  orthogonally  to,  and  tynchronously  with,  the 
motion  of  A .  B  is  delayed  by  (A'-  1)7"  t  with  respect 
to  A  to  allow  the  first  row  of  A  to  fully  enter  AODl . 

At  time  (N  -  I  ♦  k)T,  the  *th  column  of  the  matrix  B 
is  applied  to  transducers  k  through  k  +  N  of  AOD2.  In 
this  representation,  called  space  engagement,  the  ma¬ 
trix  is  folded  over  in  space  to  2N  -  1  parallel  channels 
which  require  N  time  cycles  to  be  completed.  The 
doubly  diffracted  light  is  imaged  iny  and  space  inte¬ 
grated  in  Jr  by  the  anamorphic  lens  system  following 
AOD2.  The  light  collected  on  each  photodetector 
during  each  cycle  is  the  sum  of  the  product  of  the  ele¬ 
ments  that  are  aligned  within  the  corresponding  hori¬ 
zontal  channels  of  the  two  AOOs.  As  we  will  see  in  the 
following  paragraph,  at  the  output  of  the  detector  ar¬ 
ray  we  obtain  the  product  matrix  C  *  A  X  B  in  a  time 
engagement  format. 

At  the  Nth  time  increment  the  first  row  of  matrix 
A  in  AODl  and  the  first  column  of  matrix  B  in  AOD2 
are  spatially  aligned  in  the  top  channel  of  the  system. 
The  N  local  products  alibii  are  calculated  in  parallel 
by  imaging  AODl  through  AOD2,  and  the  sum 
Cj]  *  ifij  0^6,1  is  produced  by  spatially  integrating 
all  these  products  onto  the  top  detector.  One  time  in¬ 
crement  T  later  the  second  row  of  A  has  fully  entered 
the  second  channel  of  AODl ,  and  simultaneously  the 
first  row  of  A  has  moved  one  column  away  from  the 
transducer.  In  the  orthogonal  dimension,  the  first  col¬ 
umn  of  B  has  moved  down  in  A0D2  one  channel  to  en¬ 
gage  the  second  row  of  A  and  produce  c21  *  1 

via  space  integration  onto  the  second  detector.  Concur¬ 
rently  the  second  column  of  B  has  been  applied  to 
transducers  2  through  JV  *  I ,  in  order  to  align  with  the 
first  row  of  A  and  produce  c12  *  xfijfl|,b,2  on  the 
first  detector.  In  a  similar  manner  successive  inner 
products  are  aligned,  locally  multiplied  and  globally 
accumulated  to  the  N  output  detectors.  In  a  total  of 
2N  -  1  time  increments  T ,  the  output  product  matrix 
is  produced  in  a  time  engagement  format  identical  to 
the  format  of  the  matrix  A .  Therefore  it  can  be  fed  di¬ 
rectly  back  to  the  N  transducers  of  AODl  with  no  re¬ 
formatting  or  latency.  This  allows  for  fully  efficient 
pipelining  of  iterative  matrix  algorithms,  since  after 
the  matrix  A  is  initially  loaded  into  AODl  there  are 
no  more  waiting  periods  tequired  to  load  new  ma¬ 
trices.  For  instance,  when  the  first  element  of  the  out¬ 
put  matrix  Cj  j  is  produced  the  first  row  of  A  has  been 
fully  entered  into  OADl  so  we  can  begin  entering  the 
first  row  of  the  output  matrix  in  the  top  channel  of 


AODl .  No  interference  with  subsequent  calculations 
will  occur  because  of  the  zeroes  included  in  the  space 
engagement  formatting  of  matrix  B. 

If  the  Bragg  cells  are  operated  in  the  linear  ampli¬ 
tude  diffraction  regime,  and  coherent  detection  is 
used,  then  it  is  possible  to  make  the  outputs  of  the 
photodetectors  appear  at  the  oenter  frequency  of 
AODl ,  simplifying  direct  feedback.  The  coherent  im¬ 
plementation  allows  complex  valued  matrices  to  be 
represented  by  the  magnitude  and  phase  of  the  acous¬ 
tic  pulses.  If  the  Bragg  cells  are  operated  in  the  linear 
intensity  diffraction  regime  (incoherent  implementa¬ 
tion),  then  only  real  positive  matrices  can  be  repre¬ 
sented,  but  simpler  non-interferometric  detection  can 
be  employed.  In  this  case  the  output  matrix  would 
not  be  on  a  carrier,  therefore  mixers  would  be  re¬ 
quired  to  upconvert  the  output  before  amplifying  and 
applying  to  AODl . 

3.  Frequency  multiplexed  DMAC 

An  increase  in  the  accuracy  of  this  proorssor  can  be 
incorporated  by  the  use  of  the  digital  multiplication 
by  analog  convolution  algorithm  (DMAC)  (8-12).  at 
the  expense  of  additional  complexity.  It  is  well  known 
that  the  multiplication  of  two  time  domain  waveforms 
results  in  the  convolution  of  their  Fourier  spectra. 

This  can  be  utilized  to  implement  the  DMAC  algorithm 
by  simply  multiplying  the  frequency  multiplexed  bi¬ 
nary  representations  of  two  numbers. 

The  product  z  »  x*y  of  two  M  bit  integers  x  and  y 
can  be  expressed  as  follows : 

2<Af—  1 )  M- 1  M-l  . 

z  m  S  zk2k  *(  Z/  *  2')(  y.2>) 

*«o  *  i«o  "*■ o 

2<»f-J)  ,M-l  . 

*  £  2*1  S  z,ytw),  (11 

where  xiyyt  are  the  bits  in  the  binary  representation 
of  the  integers  x  and  y  and  the  coefficients 
zk  ■  £*,.'o  *»vA_r  can  achieve  M  discrete  levels  A  time 
domain  representation  of  a  frequency  multiplexed  bi¬ 
nary  word  is  given  by  /(r)  *  Z*  .qx„  exp (jmvr).  The 
product  of  two  such  waveforms  is 
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M- 1  M- 1 

-(  E  *Jlexp(ji»wr))(  Ey„expOmwr)) 

'*■0  Am»0  / 

»  E  (  E  V/fc-JexpGAwr) 

*•0  '  »“0  7 

2(W-1) 

■  E  itexp(j*w/).  (2) 

Thus  the  weights  of  the  frequency  multiplexed  product 
waveform  correspond  to  an  A#  level  digitally  weighted 
tepresentation  of  the  product  of  the  two  binary  words. 
This  pseudo  binary  representation  can  be  channelized 
into  2M  -  1  frequency  channels,  each  centered  at  Aw. 
The  amplitude  of  each  spectral  component  can  be 
quantized  to  one  of M  levels  by  an  A/D  converter,  and 
true  binaty  representation  can  be  obtained  with  a  digi¬ 
tal  shift  and  add  register. 

Digital  multiplication  by  frequency  convolution 
can  be  incorporated  in  the  matrix  multiplier  of  fig.  1 
in  order  to  improve  the  accuracy  over  that  attainable 
with  analog  data  representation.  This  will  increase  the 
required  time  bandwidth  product  of  the  AODs  by  a 
factor  of  at  least  the  number  of  bits.  The  duration  of 
the  optical  and  acoustic  pulses  must  be  at  kast  2rr/w  s, 
to  permit  channelization  at  the  detector  output.  The 
frequency  multiplexed  binary  weighted  data  must  be 
encoded  in  phase  within  each  acoustic  pulse  so  that 
all  the  frequency  components  add  constructively. 

When  interferometric  detection  is  used,  the  RF  out¬ 
put  from  each  detector  will  be  the  coherent  sum  of  A' 
frequency  multiplexed  multilevel  binary  weighted  sig¬ 
nals,  occupying  up  to  twice  the  original  bandwidth. 
This  format  is  compatible  with  direct  feedback  to 
AOD1  without  redigitization  in  each  cycle.  Further 
iterations  would  increase  the  required  dynamic  range 
of  each  frequency  component  and  increase  the  num¬ 
ber  of  nonzero  frequencies.  By  examining  the  Fourier 
plane  of  AODI  we  can  determine  globally  the  num¬ 
ber  of  frequencies  occupied  for  matrix  A .  If  the  avail¬ 
able  number  of  frequency  bins  is  exceeded  by  p  ex¬ 
cess  Fourier  components,  then  we  can  perform  a  glob¬ 
al  pseudo  floating  point  rescaling  of  the  product  ma¬ 
trix  by  increasing  the  local  oscillator  frequency  used 
for  heterodyne  detection  by  pw.  The  lower  order  p 
bits  are  then  discarded  by  highpass  Altering  the  detec¬ 
tor  outputs,  or  with  a  Fourier  plane  aperture.  Redigi¬ 
tization  is  required  if  a  detectors  dynamic  range  is  ex- 
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ceeded  or  when  the  iterative  matrix  operation  has 
converged,  and  the  resultant  matrix  must  be  output 
in  binary  form. 

The  major  difficulty  with  incorporating  DMAC  in 
the  matrix-matrix  multiplier  is  the  complexity  of  the 
frequency  multiplexed  encoding  and  decoding,  and 
the  large  number  of  A/D  converters  required  to  con¬ 
vert  back  to  true  binary  form.  Each  of  the  N  detec¬ 
tors  must  be  channelized  into  2 M  frequency  bands, 
and  each  of  these  must  be  digitized  every  Ts  to  an  ac¬ 
curacy  of  NM.  This  could  require  as  many  as  2 NM 
A/Ds,  which  is  probably  impractical.  Alternatively, A’ 
faster  A/Ds  could  be  multiplexed  between  1M  fre¬ 
quency  channels  each  by  performing  a  conversion 
each  T/2M  s.  To  decrease  the  number  of  A/Ds  even 
further  it  may  be  possible  to  digitize  only  one  detec¬ 
tor  output  at  a  time,  and  recycle  the  matrix  for  further 
conversions  by  multiplying  by  the  identity  matrix/. 


4.  Discussion 

Acousto-optic  devices  are  attractive  transducers  for 
data  flow  optical  processors  because  of  their  sliding 
window  nature,  wide  bandwidth  (>lCHz),  large  num¬ 
ber  of  resolvable  spots  (TB>1000)  and  wide  dynamic 
range  (>60dB).  When  dealing  with  multichannel  Bragg 
cells,  however,  the  constraints  imposed  by  acoustic 
diffraction  and  electrical  crosstalk  limit  all  aspects  of 
device  performance.  Today  a  practical  limit  on  the 
number  of  channels  is  on  the  order  of  100  or  less. 
Eventually  larger  arrays  may  be  realized  through  the 
use  of  anisotropic  self  coilimation  and  effective  RF 
isolation  techniques.  The  properties  of  the  multichan¬ 
nel  acousto-optic  devices  determine  the  processing 
power  of  the  matrix  multiplier  and  to  a  lesser  degree 
the  accuracy  obtainable.  In  order  to  perform  an  A"  X  A 
analog  matrix  multiplication  this  architecture  would 
require  an  A'  channel  AOD  with  TB  *  2(2A'  -  1 ) 
(AODI ),  and  a  2A'  -  1  channel  AOD  with  TB  *  2.V 
(A0D2).  where  a  factor  of  2  was  included  for  intra 
pulse  dead  time.  If  frequency  multiplexed  binary  en¬ 
coding  is  used  with  M  bits,  then  AODI  requires  a 
TB  >M(2N  -  1).  and  AOD2  requires  a  TB  >  NM 

The  detector  array  is  composed  of  A  parallel  wide¬ 
band  photo-detectors  whose  outputs  are  bandpass  Al¬ 
tered,  combined  with  a  steering  matrix,  amplified,  and 
fed  back  to  AODI .  Single  photodetectors  can  have  a 
wide  dynamic  range  (>50dB).  but  large  monolitlucalK 
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fabricated  arrays  are  limited  by  crosstalk.  To  obtain 
the  full  dynamic  range  capabilities  of  the  detectors,  a 
powerful  optical  source  must  be  employed,  and  opti¬ 
cal  losses  minimized . 

To  obtain  an  estimate  of  the  processing  power  of 
the  space  integrating  optical  matrix  multiplier  de¬ 
scribed  in  this  letter,  let  us  consider  a  target  system 
for  multiplying  32  X  32  matrices.  AOD1  requires  32 
channels  and  a  TB  >  128,  OAD2  requires  63  channels 
and  a  TB  >  64.  If  we  assume  a  bandwidth  of  128  MHz, 
then  the  full  matrix  product  could  be  obtained  in  1  jis, 
yielding  an  analog  processing  rate  of  3.2  X  1010  mul¬ 
tiplications  per  second.  At  these  rates  the  matrix 
could  be  inverted  in  as  little  as  32  ps  with  a  direct  algo¬ 
rithm,  or  a  fraction  of  a  millisecond  with  an  iterative 
algorithm.  If  we  desire  digital  accuracy  of  8  bits  input 
and  16  bits  output,  then  the  system  parameters  be¬ 
come  much  more  stringent.  AOD1  requires  a 
TB  >  1024,  and  AOD2  requires  a  TB  >  512,  which 
could  be  accomplished  with  a  bandwidth  of  100  MHz, 
and  a  matrix  multiplication  time  of  approximately  10 
Si s.  This  yields  a  processing  rate  of  3  X  109  DMAC  mul¬ 
tiplications  per  second.  However,  to  obtain  this  addi¬ 
tional  accuracy  we  require  an  array  of  32  frequency 
channelize rs  with  16  frequency  bins  each,  and  32  8-bit 
A/D  operating  at  100  MHz,  temporally  multiplexed 
between  the  frequency  channels.  When  the  complexity 
of  the  electronic  peripherals  to  the  optical  processor 
reaches  this  high  level,  it  is  important  to  keep  in  per¬ 
spective  what  would  be  required  for  a  fully  digital  im¬ 
plementation.  For  instance,  for  the  same  numbers  as 
above  a  digital  array  of  32  8-bit  multipliers/accumula¬ 
tors  each  operating  at  1 00 MHz  can  have  the  same  pro¬ 
cessing  power  (3.2  X  109  operations/s).  In  other  words, 
comparable  electronic  hardware  is  needed  for  both  the 
digital  and  the  DMAC  optical  implementations.  The 
advantage  of  the  optical  implementation  is  in  the  re¬ 
duced  communication  rate  (320  X  106  samples/s  ver¬ 
sus  3.2  X  JO9  samples/s  in  this  example)  between  the 
high  speed  array  processor  and  the  system  buffer  mem¬ 
ory.  This  advantage  derives  from  the  2-D  parallelism 
of  optics  which  results  in  1024  parallel  multiplications 
being  performed  each  100  ns  interval. 


$.  Conclusion 

A  highly  parallel,  pipelined,  space  integrating, 
acousto-optic  processor  for  iterative  matrix-matrix 


multiplication  has  been  described.  The  architecture 
avoids  the  serial  readout  bottleneck  and  dynamic 
range  limitations  of  CCD  arrays  used  in  time  inte¬ 
grating  architectures.  The  wideband  nature  of  the  in¬ 
put  and  output  transducers  can  result  in  an  analog 
processing  rate  exceeding  30  GOPs  (Billion  multiplies 
per  second),  and  in  excess  of  3  GOPs  for  the  DM  AC 
implementation.  Additional  accuracy  can  be  incorpo¬ 
rated  by  the  use  of  the  DMAC  algorithm  and  frequen¬ 
cy  multiplexing,  but  the  improved  accuracy  is  accom¬ 
panied  by  an  increase  in  complexity  and  expense.  The 
analog  optical  processor  can  be  implemented  with  cur¬ 
rently  available  devices  and  simple  electronic  support 
circuitry.  It  provides  extremely  high  processing  power 
with  reasonably  good  accuracy  (equivalent  to  8-10 
bits)  due  to  the  high  dynamic  range  that  is  achievable 
with  non  integrating  photodiode  arrays. 
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